Origin of the giant magnetic moments of Fe impurities on and in Cs films 
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To explore the origin of the observed giant magnetic moments (~ 7hb) of Fe impurities on the 
surface and in the bulk of Cs films, we have performed the relativistic LSDA + U calculations using 
the linearized muffin-tin orbital (LMTO) band method. We have found that Fe impurities in Cs 
behave differently from those in noble metals or in Pd. Whereas the induced spin polarization of Cs 
atoms is negligible, the Fe ion itself is found to be the source of the giant magnetic moment. The 3d 
electrons of Fe in Cs are localized as the 4/ electrons in rare-earth ions so that the orbital magnetic 
moment becomes as large as the spin magnetic moment. The calculated total magnetic moment of 
M — 6.43/is, which comes mainly from Fe ion, is close to the experimentally observed value. 



PACS number: 75.20.Hr,71.70.Ej,72.15.Rn 

The problem of magnetic impurities in metals has been 
an important branch of solid state physics pL The 
Anderson impurity Hamiltonian or the Kondo impurity 
Hamiltonian has been a starting point of the model calcu- 
lations. In the side of electronic structure calculation, it 
is difficult to apply the local spin density approximation 
(LSDA) [|| in its original form, because the Coulomb cor- 
relation energy of 3d electrons is underestimated in the 
LSDA. 

Recently, giant magnetic moments are observed in the 
anomalous Hall effect measurements for the Fe and Co 
impurities on the surface and in the bulk of Cs films || . 
The deduced values of magnetic moments are ~ 7/ib and 
~ 8/iB for Fe and Co impurities, respectively. These are 
unexpectedly large magnetic moments, as compared to 
2.22/is and 1.72/1 # for pure Fe and Co metals. These 
features were attributed to an induced spin polarization 
of the neighboring Cs atoms, as in the case of 3d mag- 
netic impurities in Pd films 0. Pd is a near-magnetic 
transition element, and so Pd atoms nearby magnetic Fe 
impurity are easily spin polarized to yield giant magnetic 
moments. Whether this is also true for Fe impurities in 
Cs films is a question to be addressed in this work. 

The magnetic moments of 3c? transition metal elements 
in solids are usually determined by the spin magnetic mo- 
ments. The orbital magnetic moments are quenched by 
the crystal field and/or the hybridization effects. In con- 
trast, 4/ rare-earth elements have comparable sizes of 
the orbital magnetic moments to those of the spin mag- 
netic moments. The total magnetic moment M is given 

1/2 

by M = gj[iB [J (J + 1)] , where J is the total angular 
momentum determined by the Hund's rule: J = \L — S\ 
(L + S) for n < (21 + 1) (n > (21 + 1)), and gj is the 
Lande's ^-factor. According to the formula, a free Fe 
ion in an electron configuration of Fe 2+ will have effec- 
tively M — 6.70/is, unless the orbital magnetic moment 
is quenched [gj. This theoretical magnetic moment is 
very close to the experimentally observed M = 7\ib of 
Fe impurities in Cs jjj. Therefore, if Fe ions in Cs are 
isolated because of the large lattice constant of Cs metal, 
the main contribution to the observed magnetic moments 
will come from Fe ions themselves rather than from the 



induced spin polarization of the surrounding Cs atoms. 
In fact, previous local susceptibility data of Fe ions in 
heavy alkali-metals have been discussed in line with iso- 
lated 3d states for Fe ions The present works are 
to determine what is realized in these systems. 

We have simulated Fe impurities in Cs films using the 
linearized muffin-tin orbital (LMTO) band method ||. 
In our calculations, we assume that Fe impurities are 
substitutional and the ideal body-centered-cubic (bec) 
structure of Cs (a = 6.05 A) is maintained even after 
impurities are introduced. Although the assumptions of 
the substitutional Fe impurities and the rigid lattice are 
not valid apriori, these effects are thought to be not so 
critical ||[l0]]. The function of the basis set is adopted 
up to I = 2 orbitals for all the atomic sites. 

Since the Coulomb correlation interaction of Fe 3d elec- 
trons is expected to be significant in these systems, the 
LSDA + U method will be useful. The LSDA + U @] 
removes the deficiency of the LSDA by incorporating the 
Hubbard-like interaction term for 3d electrons. Formally, 
the LSDA + U Hamiltonian is given by 
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fc=0 

where n a mm , is the d occupation number matrix of spin er, 
F k is the Slater integral, and c (/, m : l,m') is a Gaunt 
coefficient. The conventional notations, N a = J2m n mm 
and N — N a , are used. Two main parameters in 
the LSDA + U are the Coulomb U and the exchange J 
interactions. These parameters are related to the Slater 
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FIG. 1. The LSDA DOS of lFe/5Cs ML structure. Fe 
3d-bands are extremely localized and the shape of Fe 4s 
DOS exhibit a typical 2D DOS. The spin polarization of Cs 
6s-bands is due to the hybridization with Fe bands. The spin 
magnetic moment of the first NN Cs is very small, — 0.09^ib- 
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FIG. 2. The LSDA DOS of lFe/5Pd ML structure. The 
4d DOS is shown for Pd. The hybridization of Fe 3d-bands 
with Pd Ad is large, and so the orbital magnetic moment of 
Fe is negligible. The spin magnetic moment of the first NN 
Pd is non- negligible, 0.30/j,b- 



integrals by U = F° and J = (F 2 + F 4 )/14. The ratio 
of F A /F 2 is known to be constant around 0.625 for most 
3d transition metal atoms fl2| | . We have used parameter 
values of U = 6.0 eV and J = 0.89 eV @. 

To determine the size of orbital magnetic moments, 
we have explicitly included the spin-orbit coupling in the 
Hamiltonian 14 1. In this way, our calculational scheme 
corresponds to a fully relativistic LSDA + U method. We 
have used 80 k-points sampling for the tetrahedron inte- 
gration in the irreducible Brillouin zone wedge, except 
for 40 k-points sampling in the LSDA + U calculation 
for a supercell. 

First, we have considered multilayered (ML) structures 
of Fe and Cs, which are stacks of the conventional bcc 
unit cells along the [001] direction. For lFe/3Cs struc- 
ture which contains a monolayer of Fe after three layers 
of Cs, the magnetic ordering is tested by the LSDA total 
energy calculations. The antiferromagnetic state is lower 
by only about 1 mRy than the ferromagnetic state. This 
result, however, is within the limit of the numerical pre- 
cision, and thus one cannot ascertain which state is the 



ground state. Indeed, the experiment g| identifies the 
local magnetic moment formation, but the existence of 
the long range order is not confirmed. Other experimen- 
tal tools are required to reveal the long range order in 
the system. In the following discussion, we will consider 
only the ferromagnetic ordering state. 

Figure 1 is the LSDA density of states (DOS) at each 
atomic site for lFe/5Cs ML structure. The Fe-Fe bond 
lengths are a — 6.05 A in plane and 3a = 18.15 A out 
of plane. The spin-up bands of Fe 3c? is at about —2.8 eV 
below the Fermi level Ep. In fact, it is more appropriate 
to call Fe 3d- levels rather than 3d-bands. The energy 
levels are well-defined outside of Fe 4s- and Cs 6s-bands 
so that the spin-up electrons of Fe 3d do not hybridize 
with other electrons and form local levels. On the other 
hand, the spin-down bands of Fe 3d are slightly mixed 
with Cs 6s-bands at Ep and show an energy dispersion. 
Interestingly, the shape of Fe 4s-bands is reminiscent of a 
typical 2-dimensional (2D) DOS which has a logarithmic 
divergence in the band center and becomes constant at 
the band edges. We have also found similar features for 
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Fe ions in lFe/3Cs ML. This suggests that the Fe s-s 
interaction rapidly diminishes over 2a = 12.10 A, despite 
the delocalized characteristics of s electrons, and Fe ions 
are essentially isolated in the c-direction. 

The spin splittings of Cs bands are seen for the first 
and the second nearest neighbor (NN) Cs of Fe, and it 
becomes negligible for the third NN Cs. The induced spin 
magnetic moments at Cs sites are — 0.09^_b and — 0.02^s 
for the first and the second NN, respectively. The hy- 
bridization of Cs 6s-bands with Fe bands, stronger in the 
minority band, populates the spin-down states of Cs 6s 
and drives the spin magnetic moments at Cs sites to align 
antifcrromagnetically to the spin of Fe ion. The spin and 
the orbital magnetic moments of Fe ion are, respectively, 
US = 3.46/is and /il = 1.59/is in the LSDA. When the 
LSDA + U is applied, they become /is = 3.f2//s and 
fiL = 2.87 Hb- The Coulomb correlation effects is ob- 
vious, which enlarges the 3c?-level splitting and increases 
the orbital magnetic moment. The large orbital magnetic 
moment manifests that Fe 3c? states are localized in Cs. 

A comparison with the system of Fe impurities in Pd 
will be helpful to understand Fe in Cs. In Fig. 2, we show 
the LSDA DOS for lFe/5Pd ML structure, of which the 
stacking unit is the face-centered-cubic (fee) unit cell of 
Pd. We show for Pd sites only the 4c? DOS. The exchange 
splitting of Fe 3c?-bands is still large ~ 3 eV, and the 
spin magnetic moment of Fe ion is /is = 3.12/is in Pd, 
which arc similar sizes to those of Fe in Cs. However, 
the orbital magnetic moment of Fe ion is negligible fiL = 
0.09/is in Pd. The 3d orbitals of Fe are greatly deformed 
by the hybridization with Pd, and so Fe loses most of 
ionic nature in Pd. This feature is very different from 
the case of Fe in Cs. The spin magnetic moments at 
the first and the second NN Pd sites are 0.30/is and 
0.07/is, respectively. Total sum of the spin magnetic 
moment of Fe and non-negligible spin magnetic moments 
of surrounding Pd atoms gives rise to the giant magnetic 
moment. 

For a realistic simulation of Fe impurities in Cs films, 
we have constructed a supercell, the size of which is eight 
times larger than the bec unit cell of Cs. The supercell 
has one Fe atom at the center and fifteen Cs atoms at 
the right sites. Considering that the Fe-Fe interaction 
becomes negligible over 2a — 12.10 A, the size of the su- 
percell is large enough to examine Fe impurity effects. In 
the LSDA, the properties of Fe 3c?-bands are not much 
different from those of ML structures. The spin-up bands 
of Fe 3c? shift down a little bit to higher binding energy 
(see Fig. 3). In the LSDA + U, notable differences are 
that Fe 3d states disappear near Ep by the Coulomb cor- 
relation interaction, and the energy levels of Fe 3c? are 
sharply split by the spin-orbit interaction in addition to 
the Coulomb correlation. Hence, the spin-down Fe 3c?- 
bands no longer hybridize with Cs 6s-bands. In the inset 
of Fig. 3, the angular distribution of the occupied Fe 
3c? spin-down states are plotted, based on the orbital de- 
pendent occupancies from the LSDA + U calculation. 
Although there is slight deviation from perfectly isolated 
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FIG. 3. The LSDA (upper panel) and the LSDA + U 
(lower panel) DOS's of Fe in the supercell structure. The mag- 
nitudes of the spin-up bands of Fe 3d are arbitrary. All the Fe 
bands, both of Fe 4s and 3c(, are localized. In the LSDA, one 
notices a small energy dispersion of the spin-down bands of 
Fe 3c? at Ef due to the hybridization with Cs 6s-bands. In the 
LSDA + U, Fe 3e?-levels are split by the Coulomb correlation 
and the spin-orbit interactions. Angular distribution of the 
occupied spin-down states of Fe 3d (inset) indicates that Fe 
3c? states are well-localized. 



states, the angular distribution is almost independent of 
the azimuthal angle rotation. The symmetric angular 
distribution indicates that Fe 3c? states are well-localized 
in Cs. Besides the localization of Fe 3g? electrons, the 2D 
shape of Fe 4s-bands in the ML structure is destroyed 
and Fe 4s-bands become also localized in the supercell. 

We summarize the calculated magnetic moments in Ta- 
ble E. Because of charge transfer from Cs sites, Fe ion 
gets about 0.5 (0.8) more electrons in the LSDA (LSDA 
+ U) and these electrons go to the minority bands of 
Fe 3c?. While the transferred charge reduces the spin 
magnetic moment of Fe ion, the reduced spin magnetic 
moment is compensated by the enhanced orbital mag- 
netic moment. The spin magnetic moments of Cs are 
antiparallel (parallel) to the Fe spin in the LSDA (LSDA 
+ U). The direction of the spin polarization in Cs atoms 
can be understood by the band hybridization effects. In 
the LSDA, the spin-down bands of Fe 3c? hybridize with 
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Cs 6s-bands near Ep populating the spin-down states of 
Cs 6s, whereas, in the LSDA + U, the occupied spin- 
down bands of Fe 3d are pushed down to about — 1.7 eV 
below Ep and remain as local levels. In the LSDA, the 
total magnetic moment of M = 4.62/ib is too small. It 
is because the spin-orbit coupling in the LSDA is not so 
effective to polarize orbitals. The correct orbital order- 
ing is obtained only with the strong Coulomb correlation 
interaction. The total magnetic moment of M = 6.43/ib 
in the LSDA + U is composed of Fe 94 % and Cs 6 % 
contributions [jn|. One can conclude from these results 
that the electronic states of Fe ion are almost fully local- 
ized in Cs, and the giant magnetic moment arises mainly 
from Fe ion itself rather than the induced spin polariza- 
tion of Cs atoms. In this sense, 3d electrons of Fe in Cs 
should be treated as 4/ electrons in rare-earth ions. It 
is desirable to check such a large orbital moment in this 
system with experimental tools such as the magnetic cir- 
cular dichroism (MCD) or the resonant X-ray magnetic 
scattering. The small discrepancy from the experiment 
(~ 7//b) ma Y be due to improper use of gj = 2 in the 
analysis ||, which counts only the spin contribution. 

In addition to the giant magnetic moments of Fe im- 
purities in Cs films, neither Fe coverage effects and nor 
moment change between Fe ions on the surface and in the 
bulk are observed || . In general, magnetic impurities on 
the surface of solids have larger magnetic moments than 
in the bulk. On the surface, electrons are localized due 
to lack of the degree of freedom in the surface normal 
direction. These electrons, confined on the surface plane, 
results in enhanced magnetic moments. For Fe impuri- 
ties in Cs films, Fe ions are fully localized even in the 
bulk, and thus further localization of Fe electrons does 
not occur on the surface. The negligible Fe coverage ef- 
fects can be understood considering the short interaction 
range (2a = 12.10 A) between Fe ions. The magnetic 
moment rarely depends upon the Fe concentration on the 
surface of Cs films, once excluding effects of the lattice 
relaxation and the rearrangement of surface atoms. 

In conclusion, we have found that Fe impurities in Cs 
films have a different behavior from those in Pd films, and 
that the giant magnetic moments of Fe impurities in Cs 
originate mainly from the localized Fe ions but not from 
the spin polarizations of the neighboring Cs atoms. The 

TABLE I. The calculated magnetic moments (/j,b) for the 
supercell structure. The LSDA + SO represents a relativistic 
LSDA calculation. In the LSDA + U , we have explicitly in- 
cluded the spin-orbit interaction. The magnetic moments of 
Cs are sums of all the Cs moments in the supercell. 



Method 


total M 




Fe 


Cs 


LSDA 


3.29 


Ms 


3.59 


-0.29 


LSDA + SO 


4.62 


MS 


3.54 


-0.51 






Mi 


1.60 


-0.01 


LSDA + U 


6.43 


Ms 


3.14 


0.42 






M£ 


2.90 


-0.04 



contribution from the orbital magnetic moment of Fe 3d 
electrons is comparable to that from the spin magnetic 
moment, as in the 4/ rare-earth ions. The Coulomb cor- 
relation interaction at the Fe site should be properly con- 
sidered to get consistent results with experimental data. 
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